The dynamical Casimir effect in superconducting microwave circuits 
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We theoretically investigate the dynamical Casimir effect in electrical circuits based on supercon- 
ducting microfabricated waveguides with tunable boundary conditions. We propose to implement 
a rapid modulation of the boundary conditions by tuning the applied magnetic flux through super- 
conducting quantum interference devices (SQUIDs) that are embedded in the waveguide circuits. 
We consider two circuits: (i) An open waveguide circuit that corresponds to a single mirror in free 
space, and (ii) a resonator coupled to a microfabricated waveguide, which corresponds to a single- 
sided cavity in free space. We analyze the properties of the dynamical Casimir effect in these two 
setups by calculating the generated photon-flux density, output-field correlation functions, and the 
quadrature squeezing spectra. We show that these properties of the output field exhibit signatures 
unique to the radiation due to the dynamical Casimir effect, and could therefore be used for dis- 
tinguishing the dynamical Casimir effect from other types of radiation in these circuits. We also 
discuss the similarities and differences between the dynamical Casimir effect, in the resonator setup, 
and downconversion of pump photons in parametric oscillators. 

PACS numbers: 85.25.Cp, 42.50.Lc, 84.40.Az 



I. INTRODUCTION 

Quantum field theory predicts that photons can be cre- 
ated from vacuum fluctuations when the boundary con- 
ditions of the field are time-dependent. This effect, of- 
ten called the dynamical Casimir effect, was predicted by 
G.T. Moore 1] in 1970, in the context of a cavity com- 
prised of two moving ideal mirrors. In 1976, S.A. Fulling 
et al. [2] showed that a single mirror in free space also 
generates radiation, when subjected to a nonuniform ac- 
celeration. The role of the moving mirrors in these stud- 
ies is to impose time-dependent boundary conditions on 
the electromagnetic fields. The interaction between the 
time-dependent boundary condition and the zero-point 
vacuum fluctuations can result in photon creation, for a 
sufficiently strong time-dependence [3|-|6j. 

However, it has proven to be a difficult task to ex- 
perimentally observe the dynamical Casimir effect. The 
problem lies in the difficulty in changing the boundary 
conditions, e.g., by moving physical objects, such as mas- 
sive mirrors, sufficiently fast to generate a significant 
number of photons. Although there are proposals (see, 
e.g., Ref. 0) for experimentally observing the dynami- 
cal Casimir effect using massive mirrors, no experimental 
verification of the dynamical Casimir effect has been re- 
ported to date [5|. In order to circumvent this difficulty 
a number of theoretical proposals have suggested to use 
experimental setups where the boundary conditions are 
modulated by some effective motion instead. Examples 
of such proposals include to use lasers to modulate the re- 
flectivity of thin semiconductor films [8, 9] or to modulate 
the resonance frequency of a superconducting stripline 
resonator 10], to use a SQUID to modulate the bound- 
ary condition of a superconducting waveguide [11], and 
to use laser pulses to rapidly modulate the vacuum Rabi 



frequency in cavity QED systems [12], |l3[ . 

In this paper we investigate manifestations of the 
dynamical Casimir effect in superconducting electrical 
circuits based on microfabricated (including coplanar) 
waveguides. Recent theoretical and experimental de- 
velopments in the field of superconducting electronics, 
which to a large extent is driven by research on quan- 
tum information processing |M| - [l6| . include the realiza- 
tion of strong coupling between artificial-atoms and os- 
cillators [17H19J (so called circuit QED), studies of the 
ultra-strong coupling regime in circuit QED [2fJ , single- 
artificial-atom lasing [2l[ |22| . Fock-state generation and 
state tomography [23J, l24j - Also, there has recently been 
an increased activity in studies of multimode quantum 
fields in superconducting circuits, both theoretically and 
experimentally, see e.g., Refs. [25t-t2q. and in experimen- 
tal work on frequency-tunable resonators |29h32I | . These 
studies exemplify how quantum-optics-like systems can 
be implemented in superconducting electrical circuits 
[33| , where waveguides and resonators play the roles of 
light beams and cavities, and Josephson-j unction based 
artificial atoms play the role of natural atoms in the orig- 
inal quantum optics setups. 

Here, we theoretically investigate the possibility to ex- 
ploit these recent advances to realize a system [ll| where 
the dynamical Casimir effect can be observed experimen- 
tally in an electrical circuit. We consider two circuit 
configurations, see Fig. QJc)-(d), f° r which we study the 
dynamical Casimir effect in the broadband and narrow- 
band limits, respectively. We analyze the properties of 
the radiation due to the dynamical Casimir effect in these 
systems, and we identify a number of signatures in exper- 
imentally measurable signals that could be used to distin- 
guish the radiation due to the dynamical Casimir effect 
from other types of radiation, such as thermal noise. 
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FIG. 1. (color online) Schematic illustration of the dynamical Casimir effect, in the case of a single oscillating mirror in free 
space (a) , and in the case of a cavity in free space, where the position of one of the mirrors oscillates (b) . In both cases, photons 
are generated due to the interplay between the time-dependent boundary conditions imposed by the moving mirrors, and the 
vacuum fluctuations. Here, $1 is the frequency of the oscillatory motion of the mirrors, and a is the amplitude of oscillations. 
The dynamical Casimir effect can also be studied in electrical circuits. Two possible circuit setups that correspond to the 
quantum-optics setups (a) and (b) are shown schematically in (c) and (d), respectively. In these circuits, the time-dependent 
boundary condition imposed by the SQUID corresponds to the motion of the mirrors in (a) and (b). 



The dynamical Casimir effect has also previously been 
discussed in the context of superconducting electrical 
circuits in Ref. [3J]. Another related theoretical pro- 
posal to use superconducting electrical circuits to inves- 
tigate photon creation due to nonadiabatic changes in 
the field parameters was presented in [35|], where a cir- 
cuit for simulating the Hawking radiation was proposed. 
In contrast to these works, here we exploit the demon- 
strated fast tunability of the boundary conditions for a 
one-dimensional electromagnetic field, achieved by termi- 
nating a microfabricated waveguide with a SQUID |33 |. 
(See, e.g., Ref. [36| for a review of the connections be- 
tween the dynamical Casimir effect, the Unruh effect and 
the Hawking effect). 

This paper is organized as follows: In Sec.|Hl we briefly 
review the dynamical Casimir effect. In Sec. IIIIl we pro- 
pose and analyze an electrical circuit, Fig. QJc), based 
on an open microfabricated waveguide, for realizing the 
dynamical Casimir effect, and we derive the resulting 
output field state. In Sec. HVi we propose and analyze 
an alternative circuit, Fig. QJd), featuring a waveguide 
resonator. In Sec. El we investigate various measure- 
ment setups that are realizable in electrical circuits in 
the microwave regime, and we explicitly evaluate the 
photon-flux intensities and output-field correlation func- 
tions for the two setups introduced in Sec. Ill and IV. In 
Sec. IVI1 we explore the similarities between the dynami- 
cal Casimir effect, in the resonator setup, and the closely 
related parametric oscillator with a Kerr-nonlinearity. 
Finally, a summary is given in Sec. IVIII 



II. BRIEF REVIEW OF THE DYNAMICAL 
CASIMIR EFFECT 

A. Static Casimir effect 

Two parallel perfectly conducting uncharged plates 
(ideal mirrors) in vacuum attract each other with a force 
known as the Casimir force. This is the famous static 
Casimir effect, predicted by H.B.G. Casimir in 1948 |37j |. 
and it can be interpreted as originating from vacuum 
fluctuations and due to the fact that the electromagnetic 
mode density is different inside and outside of the cavity 
formed by the two mirrors. The difference in the mode 
density results in a radiation pressure on the mirrors, 
due to vacuum fluctuations, that is larger from the out- 
side than from the inside of the cavity, thus producing 
a force that pushes the two mirrors towards each other. 
The Casimir force has been thoroughly investigated theo- 
retically, including different geometries, nonideal mirrors, 
finite temperature, and it has been demonstrated exper- 
imentally in a number of different situations (see, e.g., 
Refs. [4IH43J ). For reviews of the static Casimir effect, 
see, e.g., Refs. |USi|. 



B. Dynamical Casimir effect 

The dynamical counterpart to the static Casimir effect 
occurs when one or two of the mirrors move. The motion 
of a mirror can create electromagnetic excitations, which 
results in a reactive damping force that opposes the mo- 





Static Casimir Effect 


Dynamical Casimir Effect 


Description 


Attractive force between two conductive plates in 
vacuum. 


Photon production due to a fast modulation of 
boundary conditions. 


Theory 


Casimir (1948) 0, Lifshitz (1956) [H. 


Moore (1970) '.I], Fulling et al. (1976) [|. 


Experiment 


Sparnaay (1958) [iil, van Blokland et at. (1978) JEj, 
Lamoreaux (1997)g3], Mohideen et al. (1998) [l3j7 


— 



TABLE I. Brief summary of early work on the static and the dynamical Casimir effect. The static Casimir effect has been 
experimentally verified, but experimental verification of the dynamical Casimir effect has not yet been reported Q. 



tion of the mirror [47|. This prediction can be counter- 
intuitive at first sight, because it involves the generation 
of photons "from nothing" (vacuum) with uncharged con- 
ducting plates, and it has no classical analogue. However, 
in the quantum mechanical description of the electro- 
magnetic field, even the vacuum contains fluctuations, 
and the interaction between these fluctuations and the 
time-dependent boundary conditions can create excita- 
tions (photons) in the electromagnetic field. In this pro- 
cess, energy is taken from the driving of the boundary 
conditions to excite vacuum fluctuations to pairs of real 
photons, which propagate away from the mirror. 

The electromagnetic field in a one-dimensional cavity 
with variable length was first investigated quantum me- 
chanically by G.T. Moore [1], in 1970. In that seminal 
paper, the exact solution for the electromagnetic field in 
a one-dimensional cavity with an arbitrary cavity-wall 
trajectory was given in terms of the solution to a func- 
tional equation, known as Moore's equation. Explicit 
solutions to this equation for specific mirror motions has 
been the topic of numerous subsequent papers, includ- 
ing perturbative approaches valid in the short-time limit 
|48j , asymptotic solutions for the long-time limit |49[ , an 
exact solution for a nearly-harmonically oscillating mir- 
ror [50} . numerical approaches [5l|, and renormalization 
group calculations valid in both the short-time and long- 
time limits [52 . 15311 . Effective Hamiltonian formulations 
were reported in [54r-56j , and the interaction between the 
cavity field and a detector was studied in [53, [58| . The 
dynamical Casimir effect was also investigated in three- 
dimensional cavities 57 59], and for different types of 
boundary conditions [fjOj, |6l|. The rate of build-up of 
photons depends in general on the exact trajectory of 
the mirror, and it is also different in the one-dimensional 
and the three-dimensional case. For resonant conditions, 
i.e., where the mirror oscillates with twice the natural fre- 
quency of the cavity, the number of photons in a perfect 
cavity grows exponentially with time [62J. 

An alternative approach that focuses on the radiation 
generated by a nonstationary mirror, rather than the 
build-up of photons in a perfect cavity, was developed by 
S.A. Fulling et al. [2], in 1976. In that paper, it was shown 
that a single mirror in one-dimensional free space (vac- 
uum) subjected to a nonuniform acceleration also pro- 



duces radiation. The two cases of oscillatory motion of a 
single mirror, and a cavity with walls that oscillate in a 
synchronized manner, were studied in Refs. [63ll64| . using 
scattering analysis. The radiation from a single oscillat- 
ing mirror was also analyzed in three dimensions [651 ] . 

Table I briefly compares the static and the dynamical 
Casimir effect. See, e.g., Refs. [3|-|6| for extensive reviews 
of the dynamical Casimir effect. 



C. Photon production rate 

The rate of photon production of an oscillating ideal 
mirror in one-dimensional free space [63J , see Fig. QJa) , 
is, to first order, 
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where N is the number of photons generated during the 
time T, il is the oscillation frequency of the mirror, 
v = afl is the maximum speed of the mirror, and a is 
the amplitude of the mirror's oscillatory motion. From 
this expression it is apparent that to achieve significant 
photon production rates, the ratio v/c must not be too 
small (see, e.g., Table II). The maximum speed of the 
mirror must therefore approach the speed of light. The 
spectrum of the photons generated in this process has a 
distinct parabolic shape, between zero frequency and the 
driving frequency fi, 



— J w(0 — ui). 



(2) 



This spectral shape is a consequence of the density 
of states of electromagnetic modes in one-dimensional 
space, and the fact that photons are generated in pairs 
with frequencies that add up to the oscillation frequency 
of the boundary: u>\ + LU2 = £2> 

By introducing a second mirror in the setup, so that a 
cavity is formed, see Fig.QJb), the dynamical Casimir ra- 
diation can be resonantly enhanced. The photon produc- 
tion rate for the case when the two cavity walls oscillate 
in a synchronized manner [63l I64J , is 
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Setup 


Amplitude 
a (m) 


Frequency 

n (Hz) ' 


Photons 
n (s- 1 ) 


Mirror moved 
by hand 


1 


1 


~ lO" 18 


Mirror on a 

nano- 

mechanical 

oscillator 


10" 9 


10 9 


~10" 9 


SQUID- 

terminated 

cpw [m 


1CT 4 


10 10 


~ 10 5 



*ext(*) 



TABLE II. The photon production rates, n = N/T, for a few 
examples of single-mirror systems. The order of magnitudes of 
the photon production rates are calculated using Eq. (fl) . The 
table illustrates how small the photon production rates are 
unless both the amplitude and the frequency are large, so that 
the maximum speed of the mirror f max = aQ approaches the 
speed of light. The main advantage of the coplanar waveguide 
(CPW) setup is that the amplitude of the effective motion can 
be made much larger than for setups with massive mirrors 
that oscillate with a comparable frequency. 



where Q is the quality factor of the cavity. 

In the following sections we consider implementations 
of one-dimensional single- and two-mirror setups using 
superconducting electrical circuits. See Fig. [TJc) and 
(d), respectively. The single- mirror case is studied in 
the context of a semi-infinite waveguide in Sec. HID and 
the two-mirror case is studied in the context of a res- 
onator coupled to a waveguide in Sec. [IVJ In the follow- 
ing we consider circuits with coplanar waveguides, but 
the results also apply to circuits based on other types of 
microfabricated waveguides. 



III. THE DYNAMICAL CASIMIR EFFECT IN A 
SEMI-INFINITE COPLANAR WAVEGUIDE 



In a recent paper [11(, we proposed a semi- infinite su- 
perconducting coplanar waveguide terminated by a su- 
perconducting interference device (SQUID) as a possible 
device for observing the dynamical Casimir effect. See 
Fig. [2j The coplanar waveguide contains a semi-infinite 
one-dimensional electromagnetic field, and the SQUID 
provides a means of tuning its boundary condition. Here 
we present a detailed analysis of this system based on 
quantum network theory [6a, (67| ■ We extend our previ- 
ous work by investigating field correlations and the noise- 
power spectra of the generated dynamical Casimir radi- 
ation, and we also discuss possible measurement setups. 




FIG. 2. (color online) (a) Schematic illustration of a copla- 
nar waveguide terminated by a SQUID. The SQUID imposes 
a boundary condition in the coplanar waveguide that can be 
parametrically tuned by changing the externally applied mag- 
netic flux through the SQUID, (b) The setup in (a) is equiv- 
alent to a waveguide with tunable length, or to a mirror with 
tunable position. 



A. Quantum network analysis of the 
SQUID-terminated coplanar waveguide 

In this section we present a circuit model for the pro- 
posed device and we derive the boundary condition of the 
coplanar waveguide that is imposed by the SQUID (see 
Fig.H]). The resulting boundary condition is then used in 
the input-output formalism to solve for the output field 
state, in the two cases of static and harmonic driving 
of the SQUID. The circuit diagram for the device under 
consideration is shown in Fig. [3J and the corresponding 
circuit Lagrangian is 
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. (4) 



where Lq and Co are, respectively, the characteristic in- 
ductance and capacitance of the coplanar waveguide (per 
unit length), and Cjj and Ejj are the capacitance and 
Josephson energy of the jth junction in the SQUID loop. 
Here, $ Q is the node flux, which is related to the phase 
(p a , at the node a, as $ Q = (^o/2n)ip a , where 4>o = h/2e 
is the magnetic flux quantum. 

We have assumed that the geometric size of the 
SQUID loop is small enough such that the SQUID's 
self-inductance, L s , is negligible compared to the ki- 
netic inductance associated with the Josephson junctions 
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FIG. 3. Equivalent circuit diagram for a coplanar waveguide 
terminated by a SQUID. The coplanar waveguide has a char- 
acteristic inductance Lo and capacitance Co, per unit length, 
and it is assumed that it does not have any intrinsic dissipa- 
tion. The circuit is characterized by the dynamical fluxes $i 
and $jj. 



(&o/2ir) 2 /Ejj (i.e., a term of the form L S I 2 has been 
dropped from the Lagrangian above, where / s is the cir- 
culating current in the SQUID). Under these conditions, 
the fluxes of the Josephson junctions are related to the ex- 
ternally applied magnetic flux through the SQUID, 4> e xt, 
according to $j : i — $j,2 = ^cxt- We can therefore re- 
duce the number of fluxes used to describe the SQUID 
by introducing $,/ = ($,/.i + $j i2 )/2, and the SQUID 
effectively behaves as a single Josephson junction [69J. 

Under the additional assumption that the SQUID is 
symmetric, i.e., Cj,i = Cj t2 = Cj/2 and Ej^ = Ej t 2 = 
Ej , the Lagrangian now takes the form 

+ \c.j ($j) 2 + £./($cxt) cos U-k^- J , (5) 

with effective junction capacitance Cj and tunable 
Josephson energy 



Ej($ c 



2E.i 



$„ 
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For a discussion of the case with asymmetries in the 
SQUID parameters, see Ref. [U\ . 

So far, no assumptions have been made on the cir- 
cuit parameters that determine the characteristic energy 
scales of the circuit, and both the waveguide fluxes and 
the SQUID flux are dynamical variables. However, from 
now on we assume that the plasma frequency of the 
SQUID far exceeds other characteristic frequencies in the 
circuit (e.g., the typical frequencies of the electromag- 
netic fields in the coplanar waveguide), so that oscilla- 
tions in the phase across the SQUID have small ampli- 
tude, $j/<& <C 1, and the SQUID is operated in the 
phase regime, where -Ej^cxt) 3* (2e) 2 /2Cj. The condi- 
tion $,//"I>o <C 1 allows us to expand the cosine function 
in the SQUID Lagrangian, resulting in a quadratic La- 
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(7) 



Following the standard canonical quantization proce- 
dure, we can now transform the Lagrangian into a Hamil- 
tonian which provides the quantum mechanical descrip- 
tion of the circuit, using the Legendre transformation 
H = J2i ~§§~® i ~ £■■ We obtain the following circuit 
Hamiltonian: 



H=-Y 
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AxLq 
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and the commutation relations [$j,P,-] = ihSij and 



*i,*j 



.-, — [Pi,Pj] = 0, where Pj = -@£-. In the expres- 

sion above we have also made the identification $j = <J>i 
(see Fig. |3|). The Heisenberg equation of motion for the 
flux operator $i plays the role of a boundary condition 
for the field in the coplanar waveguide. By using the com- 
mutation relations given above, the equation of motion 
is found to be 



P 1 = C J $ 1 = -i[P 1 ,B\ 
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which in the continuum limit Ax 
boundary condition [68| 



results in the 
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0, 



(10) 



where *i(t) = *(s = 0,t), and Ej(t) = Ej [$ ext (*)] . 

This is the parametric boundary condition that can be 
tuned by the externally applied magnetic flux. Below 
we show how, under certain conditions, this boundary 
condition can be analogous to the boundary condition 
imposed by a perfect mirror at an effective length from 
the waveguide-SQUID boundary. 

In a similar manner, we can derive the equation of 
motion for the dynamical fluxes in the coplanar waveg- 
uide (away from the boundary), i.e., for $j, i > 1, which 
results in the well-known massless scalar Klein-Gordon 
equation. The general solution to this one-dimensional 
wave equation has independent components that prop- 
agate in opposite directions, and we identify these two 
components as the input and output components of the 
field in the coplanar waveguide. 



B. Quantization of the field in the waveguide 



1. Static applied magnetic flux 



Following e.g. Refs. [66|,|67|, we now introduce creation 
and annihilation operators for the flux field in the copla- 
nar waveguide, and write the field in second quantized 
form: 



+ a out ( W )e-^ a;+ ^+H.c.) 



i(—kcjX-\-{jjt) 



(ii) 



where Zq = ^Lq/Cq is the characteristic impedance. 
We have separated the left- and right-propagating sig- 
nals along the x-axis, and denoted them as "output" 
and "input" , respectively. The annihilation and cre- 
ation operators satisfy the canonical commutation rela- 
tion, [a in(ont) (uj'), al n(out) (uj")} = S(u'-w"), and the wave 

vector is defined as k u = |w|/w, where v — 1/^/CqLq is 
the propagation velocity of photons in the waveguide. 

Our goal is to characterize the output field, e.g., by 
calculating the expectation values and correlation func- 
tions of various combinations of output-field operators. 
To achieve this goal, we use the input-output formalism: 
We substitute the expression for the field into the bound- 
ary condition imposed by the SQUID, and we solve for 
the output-field operators in terms of the input-field op- 
erators. The input field is assumed to be in a known 
state, e.g., a thermal state or the vacuum state. 



C. Output field operators 

By substituting Eq. (|11[) into the boundary condition, 
Eq. (|10[) . and Fourier transforming the result, we obtain 
a boundary condition in terms of the creation and anni- 
hilation operators (for u/ > 0), 
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where 



g(u,u') 




dtEj(t)e- i(oJ -- u '^, (13) 



This equation cannot be solved easily in general, but be- 
low we consider two cases where we can solve it ana- 
lytically, i.e., when Ej(t) is (i) constant, or (ii) has an 
harmonic time-dependence. In the general case, we can 
only solve the equation numerically, see Appendix |XJ In 
Sec. |V] we compare the analytical results with such nu- 
merical calculations. 



If the applied magnetic flux is time-independent, 
Ej(t) = E°j, we obtain g(u,w') = E^J\^\S(u - u/), 
and the solution takes the form 



where 



R(w) 



a ou t(w) = H(uj) 


a in (wj, 
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Assuming that the |w| 2 Cj term is small compared to 
the other terms, i.e., that the SQUID plasma frequency 
is sufficiently large, we can neglect it in the expression 
above, and we are left with the following simplified form: 

m = _ l+^ ff R ,_ exp{2 ^ £ ff} , (16) 
Here, we have defined 



1 



cff ~ l 2n) E°jL ' 



(17) 



and assumed that k u L® s <§C 1 (this condition gives an 
upper bound on the frequencies for which this treatment 
is valid) . Figure 0] shows the dependences of Ej and L c g 
on the externally applied magnetic flux $ cx t- 

The reflection coefficient R(uj) on the simplified form 
given above [Eq. (fi"6|) ]. exactly coincides with the reflec- 
tion coefficient, — exp{2ik UJ L}, of a short-circuited copla- 
nar waveguide of length L. It is therefore natural to inter- 
pret the parameter L° ff as an effective length that gives 
the distance from the SQUID to a perfectly reflecting 
mirror (which is equivalent to a short-circuit termination 
in the context of coplanar waveguides). Alternatively, 
this can be phrased in terms of boundary conditions, 
where the mixed-type boundary condition Eq. (|10p at 
x = is then equivalent to a Dirichlet boundary condi- 
tion of an ideal mirror at x = L® s , for the frequencies 
satisfying u) -C v/L° s . See, e.g., Refs. [60, |6l| for dis- 
cussions of different types of boundary conditions in the 
context of the dynamical Casimir effect. 



2. Weak harmonic drive 

For a weak harmonic applied magnetic flux with fre- 
quency oJd, giving Ej(t) = Ej + 5Ejcos(udt), with 
SEj < E°j, we obtain 



g(u,u>') = E%I^S(uj-u J ') 

\ui\ 




(18) 
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FIG. 4. (color online) The top panel shows the normal- 
ized effective Josephson energy, -Ej^ext), and the normalized 
plasma frequency of the SQUID, cj s ($cxt), as a function of 
the external applied magnetic flux through the SQUID, $ ox t- 
The driving frequency aid, which should be much lower than 
the SQUID's plasma frequency, is also shown as a reference. 
The bottom panel shows the corresponding effective length, 
ieff(^ext)- The dashed vertical line marks the bias point used 
in the calculations, and the amplitude of the harmonic drive 
around this bias point is also indicated, by the linearized re- 
gion. The parameters used in the calculations are, if nothing 
else is specified, E°j = l.SEj, SEj = E°j/A, Ej = I c $ /(2ir), 
where I c = 1.25 fiA is the critical current of the Josephson 
junctions in the SQUID, Cj = 90 fF, v = 1.2 x 10 8 m/s, 
Z w 55 0, and ui s = 37.3 GHz, ui d = 18.6 GHz. These pa- 
rameters result in an effective length L° ff = 0.44 mm, and an 
effective- length modulation 5L g = 0.11 mm. 
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results in 



<C 1, which 



Oout(w) = R(v) a- m (u) 

S(uj,uj + ujd)e l{k " +k ^- 
S(uj,uj-ujd)e l{k " +k ^- 

where R(u) is given by Eq. (|16|) . and 

.SLeB. 



L °Ka in (uj + uj d ) + 
L °*a in (w - U) d ) + 
' )L ° a aL(u d - u), 



(20) 



S(u', u") = i-^VuTJ 1 Q(uj')Q{uj"), (21) 



where 5L eS = L° cS SEj/E j. Here, 

e = max{\S(uj,LJd — 0j)\} = 



5L e g Wd 



(22) 



is the small parameter in the perturbation calculation. 
We note that Eq. (|20| is simplified by translating it along 
the x-axis, from the point x = to the point x = L" ff 
(which is the position of the effective mirror), 

OoutM = -Oi n (w) + S((J,U + UJ d )a in (uj + UJ d ) 
+ S(u,u> -oj d )a in (uj -oj d ) 
+ S*(uj,ujd~uj)al n (uJd-uj). (23) 



It is the last term in Eq. (123)) that give rise to the 
dynamical Casimir radiation in the output field of the 
coplanar waveguide, and it appears as a consequence of 
the mixing of the aj n and a in operators due to the time- 
dependent boundary condition. Given this expression for 
ttout(w), we can in principle calculate any property of the 
output field. In Sec.|V]we discuss a number of observables 
of the output field that contain signatures of the presence 
of the dynamical Casimir part of the field described by 
the equations above. 



IV. THE DYNAMICAL CASIMIR EFFECT IN 
AN OPEN RESONATOR CIRCUIT 



Inserting this into the boundary condition above, and 
assuming that u' > 0, we obtain (after renaming u 1 — > u) 

n /in , out\ , -t t /in out\ 

= (a LJ +a iJ ) +«fe w L e ff(a w - a u ) 

ui 



2£°" 



w - uj d 



-e(u-w d )K 



UJd- w 



-&{u> d -uj)(aZ- w + ar-J 



out \ f 



UJ -\-LUd 



(a 



U!+bJ d 



h LO+UJd> 



) > 



(19) 



where 0(w) is the Heaviside step function, and where 
we also here have assumed that the SQUID is in the 
ground state. This equation cannot be solved exactly, 
because of the mixed-frequency terms (which creates an 
infinite series of sidebands around ui), but we can take 



In the previous section we discussed a setup that corre- 
sponds to a single oscillating mirror in free space. How- 
ever, we note that experimentally it might be hard to 
completely avoid all resonances in the waveguide, and 
in this section we therefore analyze the case where the 
waveguide is interrupted by a small gap at some dis- 
tance from the SQUID. Effectively this system forms an 
open coplanar waveguide resonator with time-dependent 
boundary condition, where the size of the gap deter- 
mines the coupling strength between the resonator and 
the waveguide. This system closely resembles a single- 
sided cavity in free space, with one oscillating mirror. 
A schematic representation of the circuit is shown in 
Fig. [5ja). This circuit has large parts in common with 
the circuit considered in Sec. Mil (see Figs. [5] and \S§. The 
new component is the capacitive gap that interrupts the 
coplanar waveguide, at X = 0, as shown in Fig. [SJa) and 



(b). This gap is responsible for the formation of a res- 
onator between the SQUID and the semi-infinite waveg- 
uide, at x < 0. The coupling strength between the res- 
onator and waveguide determines the quality factor of 
the resonator. This quality factor, and the correspond- 
ing decay rate, are important parameters in the following 
analysis. Note that in the limit of vanishing quality fac- 
tor this setup reduces to the setup studied in the previous 
section. 

Figure EJb) shows a lumped circuit model for the part 
of the circuit in the proximity of the capacitive gap. The 
Lagrangian for this part of the circuit is 
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(24) 



where $^ and $^ are the flux fields to the left and right of 
the capacitive gap, respectively. In the continuum limit, 
where Ax — > 0, the equations of motion for $q and $q 
result in the following boundary condition for the field in 
the coplanar waveguide on both sides of the gap: 



_L_ d$ L (a:,t) 
Lq dx 
1 9$ R (x,i) 
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(25) 



Using the field quantization from Sec. IIII[ and Fourier 
transforming the boundary condition above, results in a 
boundary condition in terms of creation and annihilation 
operators, in the frequency domain: 



~T a in( W ) 



')] 



u; 2 C c [<&(«) + aLM " «&(«) ~ aoutH] (26) 






[<M 



aJLd.; 



)] 



- oj 2 C c [<&(«) + a* t M " «&(«) - «outM] (27) 

This system of equations can be solved, and the oper- 
ators for the resonator can be written in terms of the 
operators of the semi-infinite coplanar waveguide 










,(28) 



where w c = (C c Zq)~ x is a parameter that characterizes 
the coupling strength between the resonator and the open 
coplanar waveguide. This transformation can be used to- 
gether with Eq. (|23|) , which relates the input and output 
operators of the field in direct contact with the effective 
mirror (i.e., o^(uj,x — d e g) and a^ ut (uj,x — d e s), in the 
present setup), to obtain a relation between af n (uj) and 
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FIG. 5. (color online) (a) Schematic drawing of a coplanar 
waveguide (CPW) resonator, of length d, capacitively coupled 
to an open, semi- infinite coplanar waveguide to the left, and 
terminated to ground through a SQUID to the right, (b) 
A magnification of the capacitive gap between the resonator 
and the semi-infinite waveguide, shown as a dashed box in 
(a), together with its equivalent circuit model. 



a^ ut (uj) that do not contain the resonator operators. To 
achieve this we note that the resonator operators at x = 
are related to those at x = d Q g by a simple phase factor, 
according to the transformation 



<( W ,0) 



O jfe^d of f 




a£(w,d cff ) 



(29) 



and that a^ t (w, 0) and a^(uj,0) are related to a^ ut (u) 



and a-^(w) according to the transformation in Eq. ([28 



A. Output field operators 

1. Static magnetic flux 

The resonator boundary condition on the side that is 
terminated by the SQUID is described by Eq. (|23|) . In the 
case of a static applied magnetic field the inelastic scat- 
tering by the effective mirror is absent, i.e., S(u)', to") = 0, 
and only the elastic reflections remain, 



t (u>,d e ff) = -o£(w,deff). 



(30) 



In the present setup, the SQUID is located at x = d, 
and the effective mirror is located at x = d G ff, where 
d ff = d + L® s . Thus, to write a relation between the 
input and output operator that applies on the left side 
of the resonator, we translate the boundary condition of 
the effective mirror by d c s, according to Eq. (|29l) . 

af ut (w, 0) exp{-i/c w d cf f} = -a£(w, 0) exp{ifc w d G ff}.(31) 



Since this relation applies at the point where the res- 
onator is capacitively coupled to the open coplanar 
waveguide, we can transform it using Eq. (|28j) . 
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2uJ mv ' 2u 
This equation can be rewritten as 

«outM = #res(w)a£(o>) 



where 



•RrcsM 



1 + (1 + 2iu>/u) c ) cxp{2ik LJ d c f{} 
(1 — 2iuj /lu c ) + cxp{2ifc w (i c g} 



(32) 
(33) 

(34) 



Similarly, we can apply Eq. (|28p to solve for the res- 
onator operators in terms of the input operators for the 
coplanar waveguide, 



a-outi^jdcfi) = A Tes (u)af^(u), 



where 



AresM = 



(2juj/uj c ) exp{ik u d cS } 
(1 — 2iu/uj c ) + exjp{2ik UJ d c f{ } ' 



(35) 



(36) 



The function A ros (u;) describes the resonator's response 
to an input signal from the coplanar waveguide, and it 
contains information about the mode structure of the 
resonator. From A rcs (u>) we can extract the resonance 
frequencies and the quality factors for each mode, see 
Fig. E 

The resonance frequencies, uj™ s , are approximately 
given by the transcendental equation 



tan(27< es M))= Wc /< s , 



(37) 



where ujq = 2 r Kv/d e s and n is the mode number. The 
corresponding resonance widths and quality factors are 



r, 



% Wq 



-L n 



27T V 00 n 
res 



2vr- 



' 2u> < 



(38) 
(39) 



respectively. We note that the quality factor for higher- 
order modes are rapidly decreasing as a function of the 
mode number n (see also Fig. [6]). 

Using the expressions for w^ os and r„ given above, 
the resonator response can be expanded around the res- 
onance frequencies and written in the form 



A ics {lo) 



VTn/2 



2TrT n /2-i{uj-uj T n es )' 



(40) 



and, similarly, the expression for the reflection coefficient 
of the resonator from the open coplanar waveguide is 



R TCS (uj) PS - 



r„/2 + z(q;-< cs ) 

r B /2 -»(«-«"»)' 



(41) 




U)/u>0 

FIG. 6. (color online) The absolute value of A las (cj) as a func- 
tion of the renomalized frequency co/uio, where wo = 2nv/d B fj 
is the full-wavelength resonance frequency when the resonator 
is decoupled from the open coplanar waveguide. The sequence 
of curves correspond to different coupling strengths to the 
CPW, characterized by resonator ui c values indicated by the 
labels in the figure. 



2. Weak harmonic drive 

For a time-dependent applied magnetic flux in the form 
of a weak harmonic drive, we again take a perturbative 
approach and solve the equations for a„ ut (ui) in terms of 
af n (oj) by treating S(u>',w") as a small parameter. Fol- 
lowing the approach of the previous Section, we eliminate 



the resonator variables by using Eq. 
and we obtain 



and Eq. 



t M 



i? rcs (w) 0^(w) 

5res,i(w, u + w d )a? J n (a; + to d ) 
Sns,i(w,W- w d )a? J n (a; - uj d ) 



+ S* cs . 2 (uj,uj d - u)(af n )\u) d - w), 
where R Tes (uj) is given by Eq. (|34|) . and 

S reSil (u/,u;") = S(u',u>")A ies (u')A iea (u"), 
S rcs , 2 (uj',uj") = SV,u/')4* e >')A-es(w")- 



(42) 



(43) 
(44) 



In this case the small parameter in the perturbation cal- 
culation is 



max{|5 res ,2(w,Wrf - u)\} 



8L c fi w d 1 
d e s 2 T„ 



(45) 



We note that Eq. (|42j) has the same general form as 
Eq. (|23[) . and that the only differences are the defini- 
tions of the reflection and inelastic scattering functions: 
R les (co) and S reSta {oj' , u>"). This similarity allows us to 
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(a) 



Single-mirror circuit 
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CPW 



Resonator circuit 
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■-KH 



SQUID 



FIG. 7. Schematic diagrams of the measurement setups for 
(a) the single-mirror and (b) the resonator setups in super- 
conducting microwave circuits, as discussed in Sec. Mil and 
Sec. IIVI respectively. The circulator (indicated by a circle 
with a curved arrow) separates the input and output fields 
such that only the signal from the SQUID reaches the mea- 
surement device M, and so that the input-field state is given 
by the thermal Johnson-Nyquist noise from the resistive load 
R, at temperature T. 



(a) Intensity correlation measurement setup 




(b) Quadrature squeezing measurement setup 

j r~l signal 

i 2 Q 




LO 



FIG. 8. Schematic measurement setups for (a) intensity corre- 
lations and for (b) quadrature squeezing. In (a), the box with 
t inside represents a time-delay, and in (b) the box with the 
label LO represents a local oscillator. The microwave beam 
splitter can be implemented, e.g., by a hybrid ring. The de- 
tectors are assumed to measure the intensity of the voltage 
field. In a practical experimental setup, the signals would also 
have to pass through several stages of amplification, which are 
not shown here. 



analyze both cases using the same formalism in the fol- 
lowing sections, where we calculate expectation values 
and correlation functions of physically relevant combina- 
tions of the output field operators. 



V. MEASUREMENT SETUPS 

Equation ([53) in Sec. Hill and Eq. g2j in Sec. HV| 
constitute complete theoretical descriptions of the cor- 
responding output fields, and we can apply these expres- 
sions in calculating the expectation values of any output- 
field observable or correlation function. In this Section we 
discuss possible experimental setups for measuring vari- 
ous physical properties of the output field, and we discuss 
which quantum mechanical observables and correlation 
functions these setups measure, in terms of the output- 
field operators. Below, we use Eq. (|23|) and Eq. (|42|) . and 
explicitly evaluate these physical observables for the two 
setups discussed in the previous Sections. 

The basic setups that we are considering here are 
shown schematically in Fig. [3 which also illustrates the 
concept of separating the incoming and the outgoing 
fields by means of a circulator. The input field is ter- 
minated to ground through an impedance-matched re- 
sistive load. This resistor produces a Johnson-Nyquist 
(thermal) noise that acts as the input on the SQUID. 
The circulator isolates the detector from the thermal sig- 
nal from the resistor, except for the part of the noise 
that is reflected on the SQUID. The measurement device 
is denoted by M in these circuits. 

We are interested in experimentally relevant observ- 
ables that contain signatures of the dynamical Casimir 
part of the output field, i.e., the part that is described 
by the fourth term in Eq. §Z^ and Eq. (|I5|). The most 



distinct signature of the dynamical Casimir effect is per- 
haps the correlations between individual pairs of photons, 
and such correlations could in principle be measured in 
a coincidence-count experiment. However, the physical 
quantities that can be measured in a microwave circuit 
are slightly different from those measured in the quantum 
optics regime. For instance, there are currently no single- 
photon detectors available in the microwave regime, and 
as a consequence it is not possible to directly measure 
the correlations between individual pairs of photons. In- 
stead, there are linear amplifiers 75j that can amplify 
weak signals, with very low-intensity photon-flux densi- 
ties, to larger signals that can be further processed with 
classical electronics. In addition to amplifying the signal, 
these amplifiers also add noise [7a, [77j to the signal. How- 
ever, the increased noise can often be compensated for, 
e.g., by averaging the signal over long a period of time, or 
by measuring cross-correlations in which the noise mostly 
cancel out. 

In microwave electronics, the natural fields for describ- 
ing the coplanar waveguides are the current and voltage 
fields, and these physical quantities can be readily mea- 
sured with standard equipment. Here we therefore focus 
on the voltage V(x, t) in the coplanar waveguide as our 
main physical observable. The voltage is related to the 
previously defined field creation and annihilation opera- 
tors according to 



V out (x,t) =8t$ out (a;,t) 



4tt 



ia° ut e 



-i(fc„x+a;t) + jj \ _ ( 4 g) 



The output field states described by Eq. (|23j) and 
Eq. (|42[) have voltage expectation values that are zero, 
(Vout(x,t)) = 0, but the squared voltages, i.e., as mea- 
sured by a voltage square-law detector, and various forms 
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of voltage correlations, can have non-zero expectation 
values. For example, (V ont (x,t) 2 ), (Kut(*i)K™t(*2)}, 
and (V ut(wi)V 0U t(u)2)), are all in general non-zero, and 
do contain signatures of the presence of the dynamical 
Casimir radiation. 



A. Photon-flux density 

To measure the photon-flux density requires an inten- 
sity detector, such as a photon counter that clicks each 
time a photon is absorbed by the detector. The mea- 
sured signal is proportional to the rate at which the de- 
tector absorbs photons from the field, which in turn is 
proportional to the field intensity. This detector model 
is common in quantum optics, and it is also applicable to 
intensity detectors (such as voltage square-law detectors) 
in the microwave regime, although not with single-photon 
resolution. 

The signal recorded by a quantum mechanical pho- 
ton intensity (power) detector (see, e.g., Ref. J74J]) in the 
coplanar waveguide corresponds to the observable 



I(t) ex Tr {pV^Xx^V^+Xx^t) 



(47) 



where V^ ' (x, t) are the positive and negative frequency 
component of the voltage field, respectively. In terms of 
the creation and annihilation operators for the field in the 
coplanar waveguide [see Eq. (|46j) ] , where we, for brevity, 
have taken x — 0, 



I(t) oc / dco' do/'vVu/'Tr 
Jo Jo 



pa',a u i' 



J(u'-u")t 



(48) 



and the corresponding noise-power spectrum is 

pF(-)(w)F (+) (a/) 
dio' vww' n(w,uj'), 



Sv(oj) = / du)' Tr 
Jo 

TiZp 

4^ ./„ 



(49) 



where n(oj,(j') — Tr ypa 1 (w)a(w)] . Here, Sy(oj) is related 
to the voltage auto-correlation function via a Fourier 
transform, according to the Wiener- Khinchin theorem. 
The photon-flux density in the output field, 

/>oo 

»out(w) = / doj' n out (uj,uj'), (50) 

Jo 

can be straightforwardly evaluated using Eq. (1251) . The 
resulting expression is 

"out(w) = \R(uj)\ 2 n in (uj) + \S(u,u + u d )\ 2 n in {\uj + u d \) 
+ \S(w, \ui - uj d \)\ 2 n in (\uj - Wd|) 

+ \S(LJ,LJ d -Lj)\ 2 &{LJ d -Lj) 1 (51) 

where n[ n (uj) — Tr paj n (u))ai n (ui) is the thermal pho- 
ton occupation of the input-field mode with frequency 
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FIG. 9. (color online) The output-field photon-flux density, 
n ou t{tjj), as a function of the relative mode frequency w/Wd, 
for the single-mirror setup. The solid and the dashed curves 
are for the temperatures T — 50 mK and T — 25 mK, re- 
spectively, and the dotted curve is for zero temperature. The 
red (bottom) curves show the part of the signal with ther- 
mal origin, and the blue (dark) and the green (light) curves 
also include the radiation due to the dynamical Casimir ef- 
fect. The blue curves are the analytical results, and the green 
curves are calculated numerically using the method described 
in Appendix fX] The parameters are the same as in Fig. [4] 
The presence of the dynamical Casimir radiation is clearly 
distinguishable for temperatures up to ~ 70 mK. The good 
agreement between the analytical and numerical results veri- 
fies the validity of the perturbative calculation for the param- 
eters used here. 



u, given by n™ = [exp(huj / ki,T) — 1] _1 , where T is the 
temperature and ks is the Boltzmann constant. 

The first three terms in the expression above are of 
thermal origin, and the fourth term is due to the dynam- 
ical Casimir effect. In order for the dynamical Casimir 
effect not to be negligible compared to the thermally ex- 
cited photons, we require that fc&T ^C %0Jd, where the 
driving frequency LOd here serves as a characteristic fre- 
quency for the system, since all dynamical Casimir ra- 
diation occurs below this frequency (to leading order). 
In this case it is safe to neglect the term containing the 
small factor ni n (|w +Wd|) in the expression above. Sub- 
stituting the expression for S(ui, ujd), from Eq. (I21[) . into 
Eq. (|5ip . results in the following explicit expression for 
the output field photon flux, for the single-mirror case: 



riout(w) = n- m (w) + 



5L, 



v 
SL cf{ 

v 



U \U-Ud\ Tlinflw-Wdl) 



uj(uj d - u) Q(uj d - w). (52) 
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FIG. 10. (color online) The output photon-flux density, 
Tiout(w), for the resonator setup (solid lines), as a function 
of the normalized frequency uj/ujd, for four difference reso- 
nance frequencies (marked by dashed vertical lines). Note 
that a double-peak structure appears when the driving fre- 
quency is detuned from twice the resonance frequency. For 
reference, the result for the single-mirror setup is also shown 
(red solid curve without resonances). Here, the temperature 
was chosen to be T — 25 mK, and the quality factor of the 
first resonance mode is Qo ~ 20 (cj c w 3wd), see Eq. 
The other parameters are the same as in Fig. fj] 



The output-field photon flux density, Eq. (|52p . is plot- 
ted in Fig. IH1 The blue dotted parabolic contribution to 
^out(w) is due to the fourth term in Eq. (Toll) , i.e., the 
dynamical Casimir radiation [compare Eq. ((2J)]. 

Similarly, by using Eq. (H2"j) . we calculate the output- 
field photon-flux density for the setup with a resonator, 
and the resulting expression also takes the form of 
Eq. (f5T|) . but where R(ui) and S(w',uj") are given by 
Eq. (J3J} and Eq. 03]), respectively. The resulting 
photon-flux density for the resonator setup is plotted in 
Fig. EH Here, photon generation occurs predominately in 
the resonant modes of the resonator. For significant dy- 
namical Casimir radiation to be generated it is necessary 
that the frequencies of both generated photons (u/ and 
u", where ui' + ui" — u>d) are near the resonant modes of 
the resonator. In the special case when the first resonance 
coincide with half of the driving frequency, a;™ 53 = cj^/2, 
there is a resonantly-enhanced emission from the res- 
onator, see Fig. [TO] The resonant enhancement is due to 
parametric amplification of the electric field in the res- 
onator (i.e., amplification of both thermal photons and 
photons generated from vacuum fluctuations due to the 
dynamical Casimir effect). 

Another possible resonance condition is 

W 5 CS +< S ~^- (53) 

In this case, strong emission can occur even when the 




oj/oj d 



FIG. 11. (color online) The output photon-flux density, 
n ou t{tjj), as a function of normalized frequency, for the res- 
onator setups where two modes (blue) and a single mode 
(green) are active in the dynamical Casimir radiation. The 
dashed and the dotted vertical lines mark the resonance fre- 
quencies of the few lowest modes for the two cases, respec- 
tively. The solid red curve shows the photon-flux density in 
the absence of the resonator. The blue solid curve is the 
photon-flux density for the two-mode resonance, i.e., for the 
case when the two lowest resonance frequencies add up to 
the driving frequency. The green curve shows the photon-flux 
density for the case when only a single mode in the resonator 
is active (see Fig. [10] for more examples of this case). Here, 
the temperature was chosen to be T = 10 mK, and the res- 
onator's quality factor is in both cases Qo ~ 50. The other 
parameters are the same as in Fig. [4] 



frequencies of the two generated photons are significantly 
different, since the two photons can be resonant with 



different modes of the resonator, i.e., w 
uj\ cs . See the blue curves in Fig. ITT1 



and u 



We conclude that in the case of a waveguide without 
any resonances the observation of the parabolic shape of 
the photon-flux density n ou t(w) would be a clear signature 
of the dynamical Casimir effect. The parabolic shape 
of the photon-flux density should also be distinguishable 
in the presence of a realistic thermal noise. Resonances 
in the waveguide concentrate the photon-flux density to 
the vicinity of the resonance frequencies, which can give 
a larger signal with a smaller bandwidth. One should 
note that in order to stay in the perturbative regime, 
the driving amplitude should be reduced by the quality 
factor of the resonance, compared to the case without 
any resonances. The bimodal structure of the spectrum, 
and its characteristic behavior as a function of the driving 
frequency and detuning, should be a clear indication of 
the dynamical Casimir effect. 
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B. Two-photon correlations 



The output fields described by Eq. ([23]) and Eq. (|42j) 
exhibit correlations between photons at different frequen- 
cies. This is straightforwardly demonstrated by calcu- 
lating the expectation value of the photon-annihilation 
operators at two frequencies symmetric around half the 
driving frequency, i.e., 



a ou t [ Y 



Aw) a out (— + Aw 
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V 2 
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h Aw, 

2 2 


-Aw 




"i , - ( Ud 


-Aw 



(54) 



which can be interpreted as the correlation (entangle- 
ment) between two photons that are simultaneously cre- 
ated at the frequencies Wd/2 — Aw and w^/2 + Aw, where 
Aw < LUd/2. This two-photon correlation is shown in 
Fig. [T^J for the field generated by the SQUID without 
the resonator (shown in blue), and with the resonator 
(shown in red) . Note that for thermal and vacuum states 
this expectation value vanishes for all frequencies w. This 
correlation is not directly measurable, since the operator 
combination is not Hermitian, but it serves the purpose of 
being the most basic illustration of the presence of non- 
classical two-photon correlations in the field produced 
by the dynamical Casimir effect. Below we consider two 
physically observable correlation functions that are ex- 
perimentally measurable. 



1. Second-order coherence function 

The fact that photons are predicted to be generated in 
pairs in the dynamical Casimir effect implies that the 
time-domain photon statistics exhibits photon bunch- 
ing. For instance, the measurement setup outlined in 
Fig. EJa), which measures the second-order correlation 
function 

G {2 \t) = Tr \ P V(-\0)V { -\t)V M (t)V M (Q)\ , (55) 

could be used to detect this bunching effect. For the 
output fields on the form of Eqs. ([2"3")l and (|42p . this cor- 
relation function takes the form 



G^(t) = \G^(t)\ 2 



cku w \S(uj, Ud — w)| e 



2 iur 



duj -\/w(wd — uj)R(u!d — w)5*(w, uid — w)e 4 



G^(r) 



dcjuj\S(uj,uJd — l*->)\ , 



(56) 



where R(uj) and S(lu',uj") are defined by Eqs. (|15I21[) and 
(I34I44P for the two setups, respectively. The normalized 
second-order correlation function, i.e., the second-order 
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FIG. 12. (color online) Correlations between photons at the 
frequencies w^/2 — Aw and Wd/2 + Aw, as a function of de- 
tuning Aw from half the driving frequency, for the single- 
mirror setup (blue) and the setup with a resonator (red) that 
is slightly detuned from half the driving frequency, Wd/2. The 
dashed blue curve shows the correlations for single-mirror 
setup in the presence of thermal noise, at T — 25 mK. 




FIG. 13. The normalized second-order coherence function, 
32 (t), as a function of the delay time r, for the field pro- 
duced by the single-mirror setup from Sec. IIIII (in blue), 
and for the resonator setup (in red), at zero temperature. 
The inset shows the second-order coherence function for the 
single-mirror setup without the normalization, illustrating 
that <?2 (t) > 1 and showing the oscillating behavior of gi (r) 
for large r. 



coherence function 

9 {2 \r) 



G^(i 



GW^GWir)' 



(57) 



is shown in Fig. Q2] for the single- mirror setup (blue) and 
the resonator setup (red). These coherence functions 
show clear photon bunching, since g^ 2 '(r) > 1 for a large 
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range in r. In particular, for zero time-delay, r 
coherence functions can be written as 



</ 2) (0) 



1 



0, the 



(58) 



where e is given by Eq. (|22[) for the single-mirror setup, 
and by Eq. (1431) for the resonator setup, as discussed 
in Sec. IIIII and Sec. IIV1 respectively. In both cases, e is 
small and g( 2 ' (0) 3> 1, which corresponds to large photon 
bunching. The high value of g^ 2 '(0) can be understood 
from the fact that the photons are created in pairs, so 
the probability to detect two photons simultaneously is 
basically the same as the probability to detect one pho- 
ton. For low photon intensities, this gives a very large 
second-order coherence. The decay of g^ 2 \r) is given by 
the bandwidth of the photons, which in the case without 
resonance is given by the driving frequency w<f. When 
a resonance is present, its bandwidth T determines the 
decay. Squeezed states show this type of photon bunch- 
ing [70l | . and we now proceed to calculate the squeezing 
spectrum of the radiation. 



2. Squeezing spectrum 

Another nonclassical manifestation of the pairwise 
photon correlation in the fields described by Eqs. ([25)1 
and (|4*2| is quadrature squeezing [7ll . [72J and the corre- 
sponding squeezing spectrum [73j , defined as the quadra- 
ture squeezing at a certain frequency. The quadratures 
in the frequency domain are defined by the relation 



Xg(u) 



1 



[a(w)e 



a T (w)e l 



(59) 



so that X± = Xg=Q, and X2 = Xg =7r / 2 - Experimen- 
tally, the quadratures in a continuous multimode field 
can be measured through homodyne detection, where the 
signal field is mixed with a local oscillator (LO) on a 
balanced beam splitter, resulting in a ou t(t) = (aLo(^) + 
a sig(^))/v / 2- See Fig. Efb) for a schematic representa- 
tion of this setup. The local oscillator field is assumed 
to be in a large-amplitude coherent state with frequency 
fi and phase 9, i.e., olo = M exp{— i(# + ftt)}. Probing 
the resulting output field with an intensity detector then 
provides information about the quadrature in the signal 
field, since 

I(t) = U nt (t)a out (t)) « M 2 + |q| (X e out (t)) (60) 



where 

XUt) = ~ [a sig (ty^ + 4(i) e -^+^>] . (61) 

The noise-power spectrum of the voltage intensity of the 
output field therefore gives the squeezing spectrum of the 
signal field, in the frame rotating with frequency $7: 



S 6 x (Acj) = 1+4 
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FIG. 14. (color online) The spectra of quadrature squeezing in 
the output field for a SQUID-terminated coplanar waveguide 
with (solid lines) and without (dashed lines) a resonator, as a 
function of the renormalized frequency detuning from cud/2. 
The blue (dark) and the red (light) lines correspond to the 
variances in the X e - and X e + quadrature, respectively. For 
reference, the dotted thin lines show the squeezing spectrum 
for the field produced by a parametric oscillator with a Kerr 
nonlinearity. 



where (: :) is the normally-ordered expectation value, 
and where we have normalized the squeezing spectrum 
so that S e x — 1 for unsqueezed vacuum, and S x = 
corresponds to maximum squeezing. Here, Aw is the fre- 
quency being measured after the mixing with the local 
oscillator, and it is related to the frequency u> in the sig- 
nal field as ui — fl + Acj. Hereafter, we choose il = H&. 
Evaluating the squeezing spectrum for the quadrature 



defined by the relative phase 9, i.e., X® nt (t), results in 



S x (Auj) = 1 + 2 



a ( bJd 1 a Wd a 

5 lT + Acj 'T" Aa 



2 



H 



2 



Aw 



Wd 

' 2 



Aw 



) 



+ e ™ R * (^ _ Aw) S* (^ - Aw, ^ + Aw) (63) 

where, as before, i?(w) and S(ui',uj") are defined by 
Eqs. (|15I21[) and Eqs. (J34I43I) for the two setups, respec- 
tively. 

For the single-mirror setup discussed in Sec. IIIII we 
obtain the following squeezing spectrum 



(64) 



(62) 



S x (Auj)k l-2esm(20)^/l-4( — 



where we have neglected the second term in Eq. (|63[) . 
which is one order higher in the small parameter 
S(u>', w"). Here, we can identify the relative phases 9~ = 



15 



7r/4 and 9 + = —n/4 as the maximally-squeezed quadra- 
ture (0~) and the corresponding orthogonal quadrature 
(0+). 

By using the expressions for reflection and inelastic 
scattering of the resonator setup, Eqs. (|34I43I) . in the ex- 
pression for the squeezing spectrum, Eq. (|63j). we obtain 



S e x (Aw) « 1 ± 



2e r 



l + (2^) 



(65) 



where the frequency of the first resonator mode is as- 
sumed to coincide with half the driving frequency, Wq 6S = 
cjd/2, and where we again have defined 6' ± = T 7r /4 to 
correspond to the maximally-squeezed quadrature and 
the corresponding orthogonal quadrature. The squeez- 
ing spectra for the single-mirror setup, Eq. (|64l) . and for 
the resonator setup, Eq. ([5S")l . are plotted in Fig.[T31 The 
squeezing is limited by e and e res , respectively, and it is 
therefore not possible to achieve perfect squeezing, but as 
shown in Fig.Qj] significant squeezing is still possible. In 
the single-mirror case the squeezing covers a large band- 
width, and the total squeezing (see, e.g., Ref. [73]) of the 
X g - quadrature, given by the integral of Sx(w, 9~), is 



5 Total ( 



->« «-(!-£') 



(66) 



VI. 



COMPARISON WITH A PARAMETRIC 
OSCILLATOR 



As the Q n values of the resonator considered in Sec. [TV] 
increases, its resonant modes are increasingly decoupled 
from the coplanar waveguide, and the modes become in- 
creasingly equidistant. In the limit Q n — > oo, the sys- 
tem formally reduces to the ideal case of a closed one- 
dimensional cavity (see, e.g., Refs. }lL l48l - l50l |). However, 
this limit is not realistic for the type of circuits inves- 
tigated here, because it corresponds to a regime where 
also the high-frequency modes are significantly excited, 
and this would violate our assumption that the SQUID is 
adiabatic (i.e., that the SQUID plasma frequency is the 
largest frequency in the problem). Our theoretical anal- 
ysis is also unsuitable for studying that extreme limit, 
since it implies that e rcs no longer is small. 

However, for moderate Qo _va hies, where e rcs is small 
and our analysis applies, it is still possible to make a 
comparison to a single-mode parametric oscillator (PO) 
below its threshold. The Hamiltonian for a pumped 
parametric oscillator |73J with a Kerr-nonlinearity can 
be written as 



ffi 



PO 



hu>2 t , l % 
a' a H — in 

2 2 



,—iudt. 



(at)' 



e^Va 2 



,(67) 



and where the oscillator is assumed to couple to an en- 
vironment that induces relaxation with a rate 7. The 
output field for the parametric oscillator is described by 



£» 



F{u)a™{w) + G{u)a™{-w)\ (68) 



where 



F(u) 
G(u) 



( 7 /2) 2 +q; 2 + |e| 2 
( 7 /2-zc) 2 -| e p' 

7f 

( 7 /2-*u,) 2 -|e| 2 ' 



(69) 

(70) 



see, e.g., Ref. [73j. Comparing Eq. (|68|) to the corre- 
sponding results for the dynamical Casimir effect: 



^ E H = ^ s (y+-)< CE H 



Mf- w 'f +"K CE (-") t .(7i) 



where R rcs (ui) and S ics (uj',uj") are given by Eqs. (|34l44p . 
allows us to identify relations between the parametric os- 
cillator parameters (to first order in e) and the dynamical 
Casimir parameters. We obtain 



7 = r o, 



. SL&S Ud 

' 4d e ff ' 



(72) 
(73) 



and thereby establish a one-to-one mapping between 
these systems, valid for sufficiently large Q and below 
the parametric oscillator threshold: e < 7/2, i.e., for 



5L c fi uj d 
2T 



< 1. 



-fell' 



(74) 



Using these expressions we can write a Hamiltonian that 
describes the dynamical Casimir effect in the resonator 
setup, 



Huge 



TlLU 



d 4. 

-a' a 



5L cS hud 
4d ( 



oil' 



e iui d t a 2 



e- Wdt (a^) 



(75) 



and where the capacitive coupling to the open coplanar 
waveguide induces relaxation with a rate r in the res- 
onator. This Hamiltonian picture offers an alternative 
description of the photon creation process in the dynam- 
ical Casimir effect in a resonator. This correspondence 
between the dynamical Casimir effect and a parametric 
oscillator was also discussed in e.g. Ref. |78j j . 



VII. SUMMARY AND CONCLUSIONS 

We have analyzed the dynamical Casimir radiation 
in superconducting electrical circuits based on coplanar 
waveguides with tunable boundary conditions, which are 
realized by terminating the waveguides with SQUIDs. 
We studied the case of a semi-infinite coplanar wave- 
guide, and the case of a coplanar waveguide resonator 
coupled to a semi-infinite waveguide, and we calculated 
the photon flux, the second-order coherence functions 
and the noise-power spectrum of field quadratures (i.e., 
the squeezing spectrum) for the radiation generated due 
to the dynamical Casimir effect. These quantities have 



16 





Single-mirror DCE 


Low-Q resonator DCE 


High-Q resonator DCE / PO 


Comments 


Photons created due to 

time-dependent boundary 

condition. 


The resonator slightly alters 

the mode density, compared to 

the single-mirror case. 


DCE in a high-Q resonator is 

equivalent to a PO below 

threshold. 


Classical analogue? 


No, requires vacuum 
fluctuations. 


No, requires vacuum 
fluctuations. 


Yes, vacuum and thermal 

fluctuations give similar 

results. 


Resonance condition 


no resonator 


OJ rcs = OJd/2 


OJrcs = (jJd/2 


Threshold condition 


- 


Ercs ~ Q~ 


trcs ~ Q _1 < 1 

Above threshold: nonlinearity 
dominates behavior. 


Number of DCE 
photons per second 


~ n((jj d /2)uj d 


~ n(oJrcs) r 


~ n(ajros) r 


Spectrum 
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o 
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/ \ w 


3, 
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. on resonance , 
J V w 
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o 
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* 


U! 


at T = K 




Broad peaks at resonance 

frequency uj ics and the 
complementary frequency 


Sharply 
resor 




' >. 


Broadb; 
P 


md spectrum with 
3ak at uJd/'i 


peake 
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d around the 
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TABLE III. Comparison between the dynanrical Casimir effect (DCE), in the single-mirror setup and the resonator setup, with 
a parametric oscillator (PO) with a Kerr nonlinearity. In Sec. IVII we showed that in the high-Q limit, the dynamical Casimir 
effect in the resonator setup is equivalent to a parametric oscillator. In the very-high-Q limit, the dynamics involves many 
modes of the resonator. We do not consider the latter case here. 



distinct signatures which can be used to identify the dy- 
namical Casimir radiation in experiments. 

For the single-mirror setup, we conclude that the 
photon-flux density 7i out (w) has a distinct inverted 
parabolic shape that would be a clear signature of the 
dynanrical Casimir effect. This feature in the photon-flux 
density should also be distinguishable in the presence of 
a realistic thermal noise background. 

For the resonator setup, the presence of resonances 
in the coplanar waveguide alters the mode density and 
concentrates the photon-flux density, of the dynamical 
Casimir radiation, around the resonances, which can re- 
sult in a larger signal within a smaller bandwidth. If the 
driving signal is detuned from the resonance frequency, 
the resulting photon-flux density spectrum features a bi- 
modal structure, owing to the fact that photons are cre- 
ated in pairs with frequency that add up to the driving 
frequency. The characteristic behavior of these features 
in the photon-flux density spectrum should also be a clear 
indication of the dynamical Casimir radiation. A reso- 
nance with a small quality factor could therefore make 
the experimental detection of the dynamical Casimir ef- 
fect easier. In the limit of large quality factor, however, 
the output field generated due to the dynamical Casimir 
effect becomes increasingly similar to that of a classical 



system, which makes it harder to experimentally identify 
the presence of the dynamical Casimir radiation [79( . 

For both the single-mirror setup and the resonator 
setup with low quality factor, the second-order coher- 
ence functions and the quadrature squeezing spectrum 
show signatures of the pairwise photon production and 
the closely related quadrature squeezing in the output 
field. The pairwise photon production of the dynamical 
Casimir effect has much in common with a parametrically 
driven oscillator, and in the presence of a resonance this 
correspondence can be quantified, and the two systems 
can be mapped to each other even though the systems 
have distinct physical origins. This correspondence of- 
fers an alternative formulation of the dynamical Casimir 
effect in terms of a Hamiltonian for a resonator that is 
pumped via a nonlinear medium. 
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Appendix A: Numerical calculations of output field 
expectation values in the input-output formalism 

In this section we describe the methods applied in the 
numerical calculations of the expectation values and cor- 
relation functions of the output field. Instead of taking a 
perturbative approach and solving for the output field op- 
erators in terms of the input field operators analytically, 
we can solve the linear integral equation Eq. (|12[) numer- 
ically by truncating the frequency range to [— f2, Q] and 
discretizing it in (27V + 1) steps [— wjv> -jWo = 0, ...,un], 
so that lun = fi. Here it is also convenient to define 
a(—u) = a'(cj), so that the boundary condition in the 
frequency domain reads 



= 



2tt 



cLo [e£ + a° ut ] g(u,u') 



[wfCjiaft+ag? 



vL 



fc ), (Al) 



and in the discretized frequency space takes the form 
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where we have substituted a/ —¥ uj n and uj — ¥ w m . This 
equation can be written in the matrix form 

MoutOout = M in a in => Clout = M~ ut Mi n a[ n , (A3) 

where 



M°* = - (%f) Ac 5 (c m , w „) + l -^S Un , u „ 

+ \Un\ 2 Cj5 UntU!m 

A/f in | 2lr \ a / s J\Un\ x 

M^„ = | — ) Awju ra ,w n — WreiWn 

9 J vL 



+ {UnfCjS^, 



(A4) 



(A5) 



and 



<f = (a out (^A f ),...,a out (c ),...,a out (^)) T (A6) 

a™ = (o ta («_ JV ),...,o in (wo),...,o in (« w )) T (A7) 
and, finally, where 



g(w m ,Vn) = 




2tt \ \uj. 



dtEj{t)(. 



(A8) 



which can be obtained by a Fourier transform of the drive 
signal Ej(t). 

For an harmonic drive signal we can use the fact 
that the time-dependence in the boundary condition only 
mixes frequencies that are integer multiples of the driv- 
ing frequency, and by selecting only these sideband fre- 
quencies in the frequency-domain expansion, i.e., uj n = 
u> + nuid and n = —N, ...,N, we obtain results that are 
more accurate than the perturbation results, if N > 1. 
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